Int. J. Mol. Sci. 2013, 14, 24211-24229; doi:10.3390/ijmsl41224211 



OPEN ACCESS 



International Journal of 

Molecular Sciences 

ISSN 1422-0067 

w w w . mdpi . c om/j oumal/ij ms 

Article 

Molecular Cloning and Expression of CYP9A61: 

A Chlorpyrifos-Ethyl and Lambda-Cyhalothrin-Inducible 

Cytochrome P450 cDNA from Cydia pomonella 

Xueqing Yang \ Xianchun Li ^ and Yalin Zhang 

^ Key Laboratory of Plant Protection Resources & Pest Management of Ministry of Education, 
Northwest A&F University, Yangling 712100, Shaanxi, China; E-Mail: sling233@hotmail.com 

^ Department of Entomology, University of Arizona, Tucson, AZ 85721, USA; 
E-Mail: lxc@email.arizona.edu 

* Author to whom correspondence should be addressed; E-Mail: yalinzh@nwsuaf.edu.cn; 
Tel/Fax: +86-29-8709-2190. 

Received: 8 October 2013; in revised form: 1 December 2013 / Accepted: 3 December 2013 / 
Published: 13 December 2013 



Abstract: Cytochrome P450 monooxygenases (CYPs or P450s) play paramount roles in 
detoxification of insecticides in a number of insect pests. However, little is known about 
the roles of P450s and their responses to insecticide exposure in the codling moth 
Cydia pomonella (L.), an economically important fruit pest. Here we report the 
characterization and expression analysis of the first P450 gene, designated as CYP9A61, 
from this pest. The full-length cDNA sequence of CYP9A61 is 2071 bp long and its open 
reading frame (ORF) encodes 538 amino acids. Sequence analysis shows that CYP9A61 
shares 51%-60% identity with other known CYP9s and contains the highly conserved 
substrate recognition site SRSl, SRS4 and SRS5. Quantitative real-time PCR showed that 
CYP9A61 were 67-fold higher in the fifth instar larvae than in the first instar, and more 
abundant in the silk gland and fat body than other tissues. Exposure of the 3rd instar larvae to 
12.5 mg L~^ of chlorpyrifos-ethyl for 60 h and 0.19 mg L~^ of lambda-cyhalothrin for 36 h 
resulted in 2.20- and 3.47-fold induction of CYP9A61, respectively. Exposure of the 
3rd instar larvae to these two insecticides also significantly enhanced the total P450 activity. 
The results suggested that CYP9A61 is an insecticide-detoxifying P450. 
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1. Introduction 

The codling moth Cydia pomonella (L.) (Lepidoptera: Tortricidae) is one of the most serious pests of 
pome fruits and walnuts in many countries where these crops are cultivated [1-6]. The fruit damages by 
C. pomonella include boring tunnels to the center of fruit to feed on seeds, contaminating fruits with 
their frass, and causing fruit abscission. Effective management of the codling moth mainly depends on 
chemical insecticides [2,5]. A major problem associated with chemical control of the codling moth is 
its development of resistance to insecticide groups organophosphates (OPs) [7-10] and synthetic 
pyrethroids (SPs) [11-13] in many countries. Chlorpyrifos-ethyl (OPs) and lambda-cyhalothrin (SPs) 
are the most commonly used insecticides for codling moth control worldwide and development of 
resistance to these two insecticides has been reported in the field populations of this species [3,4,10]. 
In China, however, this introduced pest has not been directly targeted by insecticides and resistance to 
insecticides has not been reported in most of the affected areas in China. 

Cytochrome P450 monooxygenases (CYPs or P450s) are a superfamily of enzymes responsible for 
detoxification of natural and synthetic xenobiotics including insecticides and plant toxins [14-19]. 
In the P450 superfamily, CYP6 and CYP9 family members are found in various insect species and play 
an important role in the detoxification of insecticides and plant toxins in insects. Exposure of insects to 
plant toxins and/or insecticides often increases the expression of a subset of P450 genes and the total 
P450 enzyme activity [15,16,20]. Involvements of CYP9 genes in insecticide detoxification and their 
inductions in response to xenobiotics have been documented in a number of insect pests including 
Helicoverpa armigera (Hiibner) [19], Aphis mellifera (L.) [17], Heliothis virescens (Fabricius) [21], 
Aedes aegypti (L.) [22], Cimex lectularius (L.) [23], and Culex quinquefasciatus (Say) [24]. 

In C pomonella, extensive studies have suggested that P450s may play a key role in detoxification of 
OP and SP due to the observed significantly higher total P450 enzyme activity in the resistant 
C. pomonella strains than in the susceptible laboratory strains [3-5,9,10]. Nevertheless, no single 
P450 gene has been isolated from this pest, not to mention the precise function of a single P450 gene in 
insecticide metabolism and resistance. To explore the potential roles of P450 genes in detoxification of 
insecticides in the codling moth, we cloned the full-length cDNA sequence of the first P450 gene, 
designed as CYP9A61, from this species. The spatiotemporal and insecticide-induced expression 
profiles of CYP9A61 revealed by real-time quantitative PCR (RT-qPCR) suggest that it is a 
xenobiotic -metabolizing P450 gene. 

2. Results 

2.1. Identification and Characterization of a Novel Cytochrome P450 Gene (CYP9A61) 

A 3' cDNA fi-agment of 992 bp and a 5' cDNA fragment of 688 bp were ampUfied by 3' and 5' RACE 
using the degenerate primers designed based on the conserved regions of the CYP9 genes (Table 1). 
A full-length cDNA of 2071 bp (GeneBank database accession number: KC832920) was then obtained 
by RT-PCR using the gene specific primers P9F and P9R (Table 1) designed based on the 
cDNA sequences of the 3' and 5' RACE fragments. The full-length cDNA sequence contains a 97 bp 
5 '-untranslated region (5'-UTR), an open reading frame (ORE) of 1638 bp, and a 357 bp 3 '-untranslated 
region (3'-UTR) (Figure 1). The ORE encodes a protein of 538 amino acids (Figure 1), with a predicted 
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molecular mass of 62.04 kDa and with a theoretical point of 7.66. The translated protein contains the 
signature motif of P450s, shares the highest amino acid sequence identity with the members of the CYP9A 
subfamily, and thus was designated as CYP9A61 by the Cytochrome P450 Nomenclature Committee. 



Table 1. PCR primers used in this study. 



Name 


Sequence (5*-3*) 


Primer used 


P9-3F1 


GGTCAAGCAGATCACCGTGaarganttyga 


3'RACE 


P9-3F2 


GGCCTTCCTGTTCTTCTTCgsnggntwyga 


3'RACE 


oligod(T) 


AAGCAGTGGTATCAACGCAGAGTAC(T)i8VN 


3'RACE 




Long: CTAATACGACTCACTATAGGGCAAGCAGTGGTATCAACGCAGAGT 


3'RACE 


10 X UPM 








Short: CTAATACGACTCACTATAGGGC 


3'RACE 


P9-5R1 


CCGCTCGGGGTCGaayttnthngg 


5'RACE 


P9-5R2 


GGCGCAGGTGGCGatnacrtcrtt 


5'RACE 


Outer Primer 


TACCGTCGTTCCACTAGTGATTT 


5'RACE 


Inner Primer 


CGCGGATCCTCCACTAGTGATTTCACTATAGG 


5'RACE 


P90RF-F 


ATGGTCTTCATTTACCTCCT 


ORE cloning 


P90RF-R 


TCAATTCTTGCGTGGCCTAAAC 


ORE cloning 


P9F 


GAAAAGTTTACATCCAACCCGAC 


Full-length cloning 


P9R 


GACGTTATAATATAACGTTTATTTCAG 


Full-length cloning 


QCYP9F 


AAATACCGCTCTGGGCATAC 


qPCR 


QCYP9R 


GATACGCCATCACCTTTACTT 


qPCR 


QActinF 


CGGCATCCACGAAACCACCT 


qPCR 


QActinR 


TGGAAGGAGCCAGTGCGG 


qPCR 



Figure 1. The full-length cDNA sequence of Cydia pomonella CYP9A61 and deduced 
amino acid sequence. The start codon, stop codon, and polyadenylation signal sequences are 
underlined. The heme-binding domain and six substrate recognition sites (SRSs) are boxed. 
Dotted lines indicate the domains on which degenerate primers were designed. 

lit DO. aa no. 

1 GmGTTTACATCCMCCCGACGTGTGITTCMCGT(nATTMTAATTTTA(lCMCATATGTTn^^ 
98 MGGTCTTCAmACCTCCTTATTGGCGCCTTCGTCGCGCTTGGCGCCTTCGTCGCGCnGG^^ 

MVFIYLLIGAFVALGAFVALGAWLYRNYTTLSHBGIKHGPFVPLLGDMAGVTFLKEHMTD 60 

278 GTCGTCAIGCGCGmMGAiUmTCCCIAACGAGAGGnGGTGGGTCGaACGAGITCCTGAaCCCT^^ 

VVMRVYESFPNERLVGRYEFLTPYVIVKDPELLKRITVKDFDHFVDRRIIVDAKAEPHF6 120 

'w^3Fi 

458 AGAAATCTCGTTTCTTTGAGAGGTGJ 




A K D L T X R Y .A .J¥..P...y...l.A TCAFGLDVHSQBEEMEFYKMSKELTHLIFETKMRFLGYRIVPHIMS 240 

P9^5R2 

|a F Iq !■ a V F e| A k| vteffsnvvlgamkqreeqhihrpdmihllmeakkgklahdtkttqddO^oo 

$RS2 SRS3 

998 GGATTTGCTACAGTAGAGGAGTtl GACGTIGGCAAAGTTAAAAACAAAACTGTGTGGTCCGACATGGATCTCGTCGCTCAGGCGGTG 

|G F A T V E E S[ D V G K V K N K T V W S D M D L V A Q A V I F L F A G Y D I V S T T M S F L L Y E L A V N P E C Q E R L 360 
SRS4 P9-3F2 
1178 CIACaAGi«anC6G6ftGaM6AMGAAaaaCfl&CGGGAA6m6&CTAIAftnCft^^ 

LQEIREHEKKHHGKFDYHSVQHMTYLDMVVSEVLRK |l P A A P Y V D R L C| V K D Y H L G K P H D K T 420 

SRSS 

1358 ACOlAAGATGnATCCTTCGTAAAGGAGAATCAGTACAAATACCGCTCTGGGCATACCAICGAGACCCCA^^ 

T K D V I L R K G E S V Q I P L W A Y H R D P K Y F P D P E K..F..D . P..E..R F S E E N K H N I N P F V Y M P |f G L G_P rT^ 480 

SRS6 Heme-binding region 

1538 ATTGGATCAAGATITGCCCTCTGCGAAGTAAAGGTGATGGCGTATCAACTACTACTCCAIATGGAGGIGACACCCTGCGAGAAGAGTATAATC 

[i~g| SRFALCEVKVMAYQLLLHMEVTPCEKSIIPITLN |m K D g| FRMAVTGGNMLRFRPRKN* 538 
1718 CATCAACATAnAAGA6TnC6IGAAGTAmACaAAAA6TmCTnAIGAAC6aTCCCTCAC(^ 

1898 GmAAGATAIKGnCGrGiaklCamiGICGaGAmGia^^ 2071 
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Amino acid sequence alignment indicates that CYP9A61 shares 58%-60% identity with the 
CYP9A members from Bombyx mandarina (Moore) (CYP9A20, 60%, ACJ05915.1; CYP9A19, 59%, 
ABQ08710.1; CYP9A21, 58%, ACJ04711.1). It also shares a high degree of amino acids identity 
with the CYP9A members from B. mori (CYP9A20, 59%, NP_00 1077079.1; CYP9A19, 58%, 
ABQ08709.1; CYP9A21, 57%, NP_001 103394.1), K armigera (CYP9A14, 59%, ABY47596.1; 
CYP9A12, 58%, ACB30273.2; CYP9A17, 56%, AAY21809.1; CYP9A18, 58%, ABB69055.1), 
and Zygaena filipendulae (L.) (CYP9A36, 57%, ACZ97417.2; CYP9A37, 54%, ACZ97407.2). 
Moreover, CYP9A61 shares more than 50% identity with CYP9A38 from Cnaphalocrocis medinalis 
(Guenee) (57%, CAZ65619.1), CYP9A39 from Spodoptera litura (Fabricius) (57%, ACV88722.1), 
CYP9A9 from Spodoptera exigua (Hiibner) (57%, ACI143222.1), CYP9A1 from H. virescens 
(55%, AAD25104.1), and CYP9A5 fromM sexta (51%, AAD51038.1). 

Figure 2. Structure-based sequence alignments of Cydia pomonella CYP9A61, 
Manduca sexta CYP9A5, Heliothis virescens CYP9A1 and Bombyx mori CYP9A20 by 
ClustalW2. The predicted a-helices and P-sheets of CYP9A61 are marked on the top of 
sequences using boxes and arrows respectively. The P-domain is red, and a-domain is blue. 
The conserved regions of P450 enzymes including FxxGxxxCxG, PxxFxPxRF, ExxR, 
WxxxR are boxed using dashed lines. 

A Helix B Helix pM pl.2 B' Helix 

■ — mmmm — — — 

cpCYF9 mvfiyixigafvaix;afvai/;awlyrnyttlshhgikhgpfvpllgdmagvtflkehmtdvvmrv^ 

MsCYP9A2 IQY IWIO'YSaYIXrrGVNYLKPLPIVGWGMMri/JKRHQAQDLDIJ.YTSFPDDRrV'GRYZMMYSKILLaZLi 

BmCY^9A20 MILLIWAWLIAAFVLF YKQAYSLFSKHGVKGFTPLPFFGNMGRIVIKMDHFSDHIQSLYDSFPZZRFVGRYZFI^PMVIIRDIZLLKKIT^ 

KvCYF9Al — MILLLTWLWI ITAVLLY FRSVYSQLSKQGVNHLPT I PVFGhnJWMVMKQZHr7DTLGRCVKAFPDDKIVGHYI»IVSPILVVLDVDTVKI^ITVKD 

C Heli x C Helix D Helix 

CpCYP9 FDHFVDRRIIVDAKAEPMFG IWLVSLRGD5;wmSiSALSPAFTSSKIRIJ^\^FMVQVGD^^IR^^ 

}isCY^9A2 SZHFIDHTSLISSSADSIFGRIFGRILFSI^GQi|WKI»13iSTLSPAFTSSiaRGMVPFMIDVNKQMVGVIKNKMQ^ 

BmCTP9A20 FZHFLDHRTIINKDTDPFFG RSLFFLaI>2rJWKDMI^^TLSPAFTSSKMKI^PFIVZVGZ^^I^^CALKQRIQZAGVGYVDIDSra 

HVCYP9A1 FZHrVDRRSFTSS-FDPIFG RGLLLIJiGD3Wig^>G;STMSPAFTSSKMRI>IVPF^ 

E Helix F Helix ^ ^xxR G Helix H Helix 



CpCTP9 AFGLDVNSQHZZ-NZFYKMSKZLTNLTFZTKMaFIX;YRIVPWIMSAFQIAVFPAian-ZFFSh^^ 

MSCYF9A2 AFGVKVNSHAZZDNZryiNGIATAIAPFRR-LISLGYSTFPSVMKKFNVKLFPTZSTTFFKNLVLRTM^ 

BmCYF9A20 AFGLKVDSITZ ZNNQFYAMGKAASTFNFRQLLIFFGIASVPiaVKILaiTLFQKilKTFFRE^ 

KvCYF9Al AFGIJCVNSQASD-HZFYVNSQAITKFKFSAFIJ<VLFFRCLPSVAQKIJ<MSLVPR£CSDYFShAA;LT^^ 

^njeli^^ ^^Jeli^^^ ^JMJeU^^^ lOielix 

CpCYP9 QDDHAGFATVZESDVGKVKNKTVWSr»mLVAQAVIFLFAGYDWSTTMSFLLYiIAWPECQZ 

}isCYP9A2 TZTDTGFAAVDZSDIGKlKVNVVWTDDDLVAQAVVFFIAGFZTISTTMTFIIJIZIJ^PDVQZiCLYAZ 

BmCY?9A20 KDSDAGFATVZZSSVGKKDINHVWTDDDLVAQAVLFFVAGFZTVSSAMTFLIJIZIJ^PZVQZKLVZZIRZNZKNNNGKFDYN 

KvCYF9Al ADADAGFATVZZSHIGaKQHNYiWTDSDLIAQAALFLFAGFDWSTSMSFU-YZIAWPDVQDRIJ/QZIRZYDZKNHGKIDYNWQS 

p2^i pl - 3 mpmp^ "'^''^ 

CpCYP9 KWPAAPVVDaiX:VKDYNLGKPNDKTTia)VILRKGZSVQIPLWAYHaDPKYFPDlPZkFDP^ 

MsCYF9A2 LWPPLVSLDaI/:IIa)YNLGRANRKCKKGFIIRKGZSVAVPVYSIHHDPZYFPziPYRFNPDRSSZG^^a^DIKE^YMP;^ 

BmCYF9A20 LWPPVIALDHMCVKDYNLGKPNDKSKZDFIiaKDVAVGIPWGIJiRDPZFFPb?PLKFDPZRF5ZZNKHNIKPFSYMPlFGLGPRNCIGisRF 

HVCYP9A1 LWPPAAVVDRVCVKDYNIGRPNKKATKDLIIHTGQAVAISPWLFHRNPKFFPZjPj^^ 

P3 P4-1 P4-2 PxxFxPxRF FxxGXXXCxG 

CpCYF9 QLLLHMZVTPCZKSIIPITLNMKDGFRMAVTGGNWLRFRPRKN 

MSCYF9A2 ZVLQHFZVKVCZKTTIPARLDTTT-LSMGIDGGHWVRFKLRH- 

BaiCYP9A20 QLLQHMEISPCZKTCIPSKLSKET-FNLRLZGGHWIRLKIRN- 

HvCTP9Al LL 



The protein contains 16 a-helices, 9 P-sheets and the loops are linked (Figure 2). The second structure 
of this enzyme is composed of a P-domain (denoted in red: A helix, B helix, pi, P2), and an a-domain 
(denoted in blue). Motif prediction by scanning the Prosite database shows the CYP9A61 contains a 
number of typical characteristic features of cytochrome P450 (Figure 2). For instance cytochrome 
P450 heme-binding domain FxxGxxxCxG (residues 473-482, in which C is the heme iron ligand) is at 
the C-terminal that is conserved in most of insect P450 enzymes. The typical aromatic motif PxxFxPxRF 
(meander) on the K-K' loop is found at the position of 449-459 aa in this protein. The ERR-triad formed 
by the ExxR (residues 392-395) motif and PxRF (residues 454-457) motif may form a set of salt-bridge 
interactions to assure the cysteine pocket is at the right position to make sure the heme binds to 
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protein [25,26]. The conserved WxxxR (131-135) region located near the A^-terminal in the C Helix may 
be involved in eliminating the charge of propionate side chains of the heme [27]. Six substrate 
recognition sites (SRSs) were also determined (Figure 1) based on the CYP9A61 secondary structure 
elements [26]. SRSl is located in the loop B'-C helixes, SRS2 and SRS3 are situated in the helices F and 
G connecting the loop. In general, SRS4 lies in the center portion of helix I; however, it is located in the 
loop H-I helixes in this protein. SRSS is found at the A^-terminal of P strand 1-4, and P strands 4-1 and 4-2 
protrude into SRS6. The CYP9A61 was predicted to have no signal peptide in the protein A^-terminal but 
contains a transmembrane domain (residues 2-24) at the A/-terminus (Figure 1). 

2.2. Phylogenetic Analysis 

To examine the relationship of CYP9A61 with members of the CYP9 family from Lepidoptera 
deposited in the GenBank database, a phylogenetic tree was constructed using the neighbor-joining 
method (Figure 3). The analysis shows the CYP9A61 is located on the same branch with 
Danaus plexippus (L.) CYP9 andM sexta CYP9A5 with 48% bootstrap support, indicating CYP9A61 
shared closest relationship withZ). plexippus CYP9 midM. sexta CYP9A5. Although the Bombycidae's 
CYP9s share the highest amino acid identity with CYP9A61 (57%-60%), they are located in another 
branch with the CYP9s from Noctuidae, Pyralidae and Zygaenidae. In this tree, the transcripts of some 
CYP9 members were observed to be induced by insecticides. For instance, deltamethrin increased the 
expression levels of C7/^P^77 and CYP9A12 fxomH. armigera [26]. The mRNA level of CYP9A1 from 
H. virescens was increased through exposure to thiodicarb [21]. 

Figure 3. Phylogenetic relationship of C. pomonella CYP9A61 with 22 CYP9s from 
Lepidoptera. CYP9A61 (GenBank accession number: KC832920) is boxed. This un-rooted 
phylogenetic tree was constructed using the neighbor-joining method. Nodes indicate 
bootstrap calculated with 1000 replications support. 




^^B.mori CYP9A19 ABQ18318.1 



66 



43 



l_ B.mori CYP9A 1 9 ABQQH 709. 1 

r B.mandahnagi CYP9A20V1 ACJ05915.1 

^^B.moh CYP9A20NP001 077079.1 
Cmedinalis CYP9A14 CAZ65619J 



Zfilipendulae CYP9A37 ACZ97417.2 



82 



86 
61 



- .V. exigua CYP9A 9 BA G7 14 10. 1 
M.brassicae CYP9 AAR265 1 8. 1 
'Karndgera CYP9A17 EU54 1247.1 



93 



m^Karmgera CYP9AI2 EU54 1248.2 
— B.mori CYP9A22 ABQ08707.1 



100 




Zfilipendulae CYP9A36 ACZ974 1 7.2 



jRzea CYP9A14 ABH09253.1 



'\^H.armgera CYP9A14 EU327674.1 
Kvirescens CYP9A1 AAC25787.1 



99 



Pxylostella CYP9G4 ACH88357.1 

- P.xylosteIIa CYP9G2 BA 071413.1 




\Cydia pomonella KC832920 \ 
^D.plexippus CYP9 EHJ65261.1 
M.sexta CYP9A4AAD5 1038.1 
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2.3. mRNA Expression of CYP9A in Response to Chlorpyrifos-Ethyl and Lambda-Cyhalothrin Exposure 

In this study, we investigated if chlorpyrifos-ethyl and lambda-cyhalothrin are inducers and if 
CYP9A61 is involved in detoxification of insecticides. Third instar larvae were treated with 
chlorpyrifos-ethyl and lambda-cyhalothrin, and their mRNA expression levels were detected by 
RT-qPCR. After treatment with chlorpyrifos-ethyl on the dorsum of the larvae, no significant difference 
in CYP9A61 expression was observed in the first 48 h. The 60 h chlorpyrifos-ethyl exposure led to a 
2.2-fold up-regulation of CYP9A61 at a concentration of 12.5 ng when compared with the 
control (Figure 4A). The CYP9A61 transcript level was strongly induced after exposure to 
lambda-cyhalothrin as time progressed from 12 to 36 h (approximately 1.32-3.47 fold), and thereafter 
dropped to 2.60 fold after 60 h of lambda-cyhalothrin exposure (Figure 4B). 

Figure 4. Induction of C. pomonella CYP9A61. RT-qPCR analysis relative expression of 
CYP9A61 treated with 12.5 ng chlorpyrifos-ethyl (A) and 0.19 ng lambda-cyhalothrin (B). 
Expression level was first normalized to the reference gene hetsi-actin. The normalized value 
of each sample was then divided by a specified control (the value of CYP9A61 in acetone 
treated control), and the ratio was applied to relative expression analysis. The results are 
shown as the mean ± SE. The error bars show the ranges of standard errors. Letters on the 
error bars indicate the significant differences by ANOVA analysis (p < 0.05). 



3n 



2 2- 



Chlorpv rifos-ethyl 



■ 
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X 



-r 

Control 12h 24h 36h 48h 60h 
Exposure time 



4n 



3- 



2- 



1- 
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b b 



-T 



T- 



Control 12h 24h 36h 48h 60h 
Exposure time 



2.4. Developmental Expression Profiles 

To determine the relative expression patterns of CYP9A61 mRNAs at the various developmental 
stages, a RT-qPCR was conducted. The CYP9A61 transcript was detectable in all developmental stages, 
and was significantly higher in the fifth instar than at other ages (Figure 5 A). The results showed that 
no significant difference in CYP9A61 expression was observed from first to third instar. There was a 
rapid increase of CYP9A61 expression at the fourth instar of C. pomonella, and it reached a peak in 
fifth instar larvae (67-fold of the first instar). The gene expression level dropped at the pupal stage 
(10.85 and 6. 15 -fold of the first instar for male and female, respectively) and finally increased in adults 
(44.44 and 26.94-fold in male and female, respectively). Interestingly, the expression of CYP9A61 was 
not different between the male and female codling moth in pupal and adult stages. 
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Figure 5. Relative expression level of C. pomonella CYP9A61 in different developmental 
stages (A) and tissues (B). Expression level was normalized using (3-actin as the standard. 
The normalized value was applied to relative expression analysis. The results are shown as 
the mean =t SE. The error bars show the ranges of standard errors. Letters or asterisks on the 
error bars indicate the significant differences by ANOVA analysis {p < 0.05) or by Student's 
Mest 0.001). 
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2.5. Tissue Specific Expression Profiles 

The mRNA expression levels of the CYP9A61 in the head, cuticle, silk gland, midgut, fat body, and 
Malpighian tubule dissected from third instar larvae were investigated (Figure 5B). The results showed 
that CYP9A61 was expressed in all tissues (Figure 5B). Figure 5B shows that the CYP9A61 has a high 
level of mRNA expression in silk glands (127.78-fold higher than in the cuticle) and fat body 
(384.00-fold higher than in the cuticle), while the gene has relative low expression in midgut, head and 
the Malpighian tubule, at 9.87-, 3. 61-, and 12.58-fold of the cuticle (Figure 5B). 

2. 6. Enzyme Activity 

In order to determine whether chlorpyrifos-ethyl and lambda-cyhalothrin could influence the 
mixed- function oxidase (MFO) activity, the MFO activities using /^-nitroanisole as substrate are shown 
in Figure 6. The between- subjects effects (insecticide exposure and exposure time) were analyzed using 
a factorial ANOVA (Table 2). A significant difference in j^-nitroanisole demethylase (PNOD) activity 
was observed after treatment with chlorpyrifos-ethyl (Table 2, Figure 6). The PNOD activity was no 
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different between the treated group and control after 12 h treatment, but significantly increased after 
24 and 36 h treatment, and then fell after exposure for 48 h. The PNOD mean activity for 24 h 
chlorpyrifos-ethyl-exposure is 11.55 pmole min~^ mg protein"^ and the value is 7.14 pmole min~^ mg 
protein^^ for the control, which constitutes a 1.62-fold increase. A weak activity (0.83, and 0.83-fold 
compared to control) was observed for the experimental group at 48 h and 60 h, with the value of 
6.75 and 6.66 pmole min~^ mg protein~\ respectively. Higher activity was detected for the longer time 
inductions (72 and 96 h), with 1.14 and 1.20-fold of PNOD activity increase when compared with 
the control. 

Figure 6. MFO activities of C. pomonella treated with 2.5 ng chlorpyrifos-ethyl and 
0.19 ng lambda-cyhalothrin using PNOD as substrate. The results are expressed as pmole 
/^-nitrophenol min~^ mg of protein~\ The error bars represent the standard error across 
three replicates. Asterisks on the error bars indicate the significant differences between 
chlorpyrifos-ethyl or lambda-cyhalothrin treatment and control by Student's ^test. 
p< 0.05; -^-^ p< 0.01; p < 0.001. 
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Table 2. Variance analysis of factors influencing PNOD activity induced by two insecticides 
on Cydia pomonella larvae (between- subjects effects). 



Source 




Mean square 


F 


p" 


Chlorpyrifos 










Insecticide exposure 


1 


10.694 


22.775 


0.000 


Exposure time 


7 


44.106 


99.932 


0.000 


Insecticide exposure x exposure time 


7 


9.781 


20.831 


0.000 


Error 


80 


0.470 






Lambda-cyhalothrin 










Insecticide exposure 


1 


288.956 


745.795 


0.000 


Exposure time 


7 


138.226 


150.066 


0.000 


Insecticide exposure x exposure time 


7 


44.372 


47.087 


0.000 


Error 


80 


0.921 







^ Degree of freedom; Distinctness index,/? < 0.05 indicates that there are significantly different. 
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For the lambda-cyhalothrin exposed larvae, the main interaction effects for insecticide exposure and 
exposure time and insecticide exposure and exposure time were significant (Table 2). Significantly 
higher PNOD activities were exhibited after exposure to lambda-cyhalothrin in the previous 72 h 
compared with the control, and then declined to the same level as the control (Figure 6). The results 
indicate that the PNOD activity was 2.09-, 2.05-, 1.98-, 1.81-, 1.90- and 1.31 -fold higher than the control 
after exposure to lambda-cyhalothrin as time progressed from 12 to 96 h, respectively. No difference 
was observed after 96 h of lambda-cyhalothrin exposure. 

3. Discussion 

The great interest in insect P450s is due to their role in detoxification of natural and synthetic 
endogenous and exogenous compounds including insecticides and plant toxins [14-19]. CYP9A61 
cloned from C pomonella larvae in the present study represents the first P450 gene isolated from 
this species. Sequence alignment demonstrates that CYP9A61 shares 57% and 56% homology with 
K armigera CYP9A12 and CYP9A17, respectively. In addition, the SRSl, SRS4 and SRS5 are highly 
conserved with the CYP9 members above (Figure 7). In H. armigera, CYP9A12 and CYP9A17 were 
potentially involved in xenobiotic metabolism [26]. It is noteworthy that CYP9A61 shares a 55% amino 
acid identity with H. virescens CYP9A1 and 51% with CYP9A5 firomM sexta\ these two P450s were 
found to be induced by thiodicarb [21] and clofibrate [28], respectively. Therefore, we infer that the 
CYP9A61 in codling moth could also be induced by insecticide exposure. Based on the high sequence 
identity and second structure similarity, we infer that the CYP9A61 have the same function in metabolism 
of endogenous and exogenous compounds including insecticides, fatty acids and plant toxins. 

We investigated the CYP9A61 transcription levels in different development stages. Interestingly, 
the results show the highest value in fifth instar larvae and levels then dropped in the pupa and adult 
which is congruent with previous studies in a NADH-cytochrome b5 reductase from H. armigera [29], 
suggesting that the CYP9A61 expression profile during the life cycle is consistent with that of 
co-regulated functional elements of the microsomal drug oxidation system. Furthermore, the expression 
pattern is also in line with that in the Diptera insect, C. quinquefasciatus [24]. No significant difference 
in expression of CYP9A61 was observed among males and females in both pupa and adult stages, 
suggesting that CYP9A61 expression may not be related to gender in these stages. However, CYP9A12 
from H. armigera was found to have a higher expression level in both the pupal and adult male than in 
the female [26]. Furthermore, CYP transcription profiles of some CYP6s from Dipteran species, such as 
Drosophila melanogaster (Meigen) [30] and CYP6Z1 irom Anopheles gambiae [31] where there is a 
significant difference of gene expression among males and females, may be related to the gender. 
Unfortunately, sexing the larvae is unfeasible and whether CYP9A61 is related to gender in our larvae 
was uninvestigated. The gene expression patterns of P450s in different developmental stages were 
related to the biological and physiological function of P450s. Our results reveal that CYP9A61 
expression levels in the whole body of the codling moth are higher in the active feeding stages 
(fourth and fifth instar larvae) than non-feeding stages (pupa and adult). This is somewhat consistent 
with CYP4 genes from fat body and midgut of M sexta [32]. In that case, the expression levels of 
CYP4 genes were higher in active feeding, midwandering, prepupal, and pupal stages. This expression 
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pattern suggests that CYP9A61 is potentially responsible for detoxification of xenobiotic compounds 
including insecticides or plant toxins suffered during feeding. 

Figure 7. Comparison of C pomonella CYP9A61 and CYP9A12 (EU541248.2) and 
CYP9A17 (EU541247.1) from Helicoverpa armigera, SRSl, SRS4 and SRS5 are boxed. 
Asterisks (*) indicate identical residues among three sequences; colons (:) indicate residues 
with conserved substitutions; dots (.) represent residues with weakly conserved residues in 
the three sequence alignment. 



BACmAl 7 MILALVHVKVLIAVAVL yiAQVYSRmYGVKQlltFVPILGMMmLIJC^ 

BMCyP9A12 MILVLVNVAVLITVSVL YIAQIYSRmireVKHFRFVPIXGmflBMVUOQDHLI 

CpCyP9A€l MVFIYIXIGAFVALSAFVALGmyiOIYTTLSHHGIIffl 

CoDsensus : :* : *; :*::*:*: • • ;;*; ; •* ***** **•**••** 



HaCYP9A17 
SaCyP9A12 
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CoDseDSUS 



SRS 1 



IKELWIRDIELAKKIAVKDFEHFLDHRSMFSS-SDSFFS RNLFSLK SQEWm^RSTLSPAFISSKMRMMVPFMVEVGDQ 
ISEMWIRDLELVKKIAVKDFEHFLDHRSIFSS-SDSFFS RNLFSLK 3QEWKDMRSTLSPAFTSSKMRMKVPEHVEVGDQ 
LTPYVIVKDPELLKRITVKDFDHFVDRRI IVDAKAEPMFg |r!ILV5LBP DEWKDMRSALSPAFTSSKIRLMVPFMV^ 
:. *:::* ** *:*•****•**•*•* •.: :••*.*****:*•*******•*********:*•******•**** 



BaCYPSAl 7 MAAIKNKIKESGNGYIDICCKDLmYANDVIASCAroLKVDSHHETI^ 

HaCyP9A12 lftlD2^KKKIKES6NGyiDLECia)LTTRyANimASCAF6LK\a>SHNET 

CpCYP9A€l MIRALKTKIEQSQVDYIDLEJUa)LTTRYANDVIATa^LDVNSQHE-ENEFY!^ 
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HaCYP9A17 
SaafP9A12 
CpCyP9A61 
Consensus 

BaCYP9A17 
HaCYP9A12 
CpCyP9A€l 
Consensus 

BaCyP9A17 
BaCyP9A12 
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SRS 4 



KILKLOFLSCAAKKFFRNLVLDTMKNREXNHIIRPDMIHLIifEAKKGKLTHEEI lAGFAIVEES IIVGQKEITRVW 
KIUaJ>ELSEAAKKFnaiLVLDIMKNRELI«IIRPIXaaLU«^^ lAGFATVEES IVVGQKEITRVH 

SAFQLAVFPJtfgVIErFSNVVLSMflCggEO^ 
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TDEDLIAQAVLFFIAGFETVSSOlSFLLyELAWPDVQERLAQEIKENDAKNGGKFDFNSIQNLQYMD^^ 
TDEDLIAQAVLFFIAGFETVSSCaiSFLLYELAVNPDVQERLAQEIKENDAKNGGKFDFNSIQNLQYMD^^ 
SDMDLVAQAVIFLFAGYDTVSTTMSFLLYELAVNPECQERLLQEIRENEKKNNGKFDYNSVQ^ 

•* **•****•*••**••***• ************* **** ***;*** **,****•**•*•• *********** **. 

SRS 5 

^KDYNLGKPHIHC2aCH0FIVI«STGISIPim^FHEU)PQFm 
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HaCYP9Al 7 CIGSRFALCEMKVMAYQILQHMEVSPCERTCIPAiaDTET-FlOlRLKGGHWLRFRPRQ- 

HaCYP9A12 CIGSRFALCEMKVMAYQILQHMEVSPCERTCIPAKLOIET-FNMRLKGGHNLRFRPRQ- 

CpCyP9A€l CIGSREUXXViCVMAYQLLLmiEVTPCEKSIIPITUDOa)GF^ 

Consensus **********>******■* ****:***•> ** .*• • *.* ■.**•*******• 



To better understand the function of CYP9A61, the distribution of CYP9A61 in the tissues of the 
codling moth was determined. It has been reported that tissue-specific expression of detoxification 
related genes might illustrate some of the functions of their biological and physiological roles [31]. 
We found similar developmental expression profiles of CYP9A61 with NADH-cytochrome b5 reductase 
from H. armigera, but the tissue-specific expression pattern of CYP9A61 is inconsistent with the 
H. armigera NADH-cytochromeytochroreductase [29]. To our knowledge, the fat body in insects is a 
key tissue associated with energy metabolism, detoxification of xenobiotics and intermediaries [33,34]. 
In our study, the high expression level observed in fat bodies further indicates that CYP9A61 may be 
associated with detoxification of xenobiotic compounds. Silk glands, an organ which produce threads of 
silky material, where CYP9A61 is abundantly expressed, indicates that CYP9A61 has a potential role in 
protection of codling moth from mechanical injuries. However, to determine the accurate function of 
CYP9A61 in detoxification of xenobiotic compounds, the CYP9A61 expression level of specific-tissues 
of codling moth exposure to insecticides should be further investigated. 

In insects, substantial experimental effort is being focused on the P450s in response to 
plant allelochemicals and other secondary metabolites to elaborate plant-insect co-evolution 
interactions [18,35]. The specific P450 genes acting in response to insecticide exposure are less 
known [36], especially in the codling moth. A previous study demonstrated that the MFO activities for 
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all stages of whole larvae in a pyrethroid resistant codling moth population were significantly increased 
compared to a susceptible strain [37]. However, inducibility is one of the characteristics of P450s; many 
of them can be induced by xenobiotic compounds [15,38]. In a thiodicarb selected H. virescens strain, 
the expression of CYP9A1 from whole body was 29- and 13 -fold higher than in two susceptible strains [21]. 
Elevated CYP9 mRNA transcripts for all developmental stages were observed in pesticide-exposed 
C lectularius collected under repeated insecticide treatments in an apartment in Columbus, OH, 
USA [23]. In A. mellifera, CYP9Q3 in the midgut of 3-d-old worker bee was slightly induced (1 .5-fold) 
by tau-fluvalinate [17]. A similar phenomenon has been reported for CYP9M8 (2.5-fold) 
and CYP9M9 (2-fold) in permethrin-exposed fourth instar A. aegypti larvae [22]. Furthermore, in 
H. armigera, deltamethrin induced CYP9A12 and CYP9A17 expression in the midgut whereas it 
repressed both CYP9A12 and CYP9A17 transcripts in fat body [26]. Based on this evidence, we 
hypothesize that lambda-cyhalothrin and chlorpyrifos-ethyl can also induce CYP9A61 in the codling moth. 

In this study, the response of CYP9A61 to chlorpyrifos-ethyl and lambda-cyhalothrin exposure was 
investigated. The CYP9A61 transcript is significantly induced by lambda-cyhalothrin (Figure 4B) 
using the dosage selected, and the same tendency is shown in the enzyme activity assay (Figure 6). 
Many CYP genes belonging to the CYP9 family were induced by synthetic pyrethroid insecticides. 
In D. melanogaster, xenobiotic inducibility of the P450 genes is associated with insecticide 
resistance [30]. In pyrethroid-selected A. aegypti strains, more than 13 CYP9s were observed 
up-regulated [39]. Gonzalez et al, observed a high level of CYP in fenitrothion- and pyrethroid-resistant 
populations [40]. Yang et aL demonstrated that CYP9A12 and CYP9A14 from H. armigera are 
up-regulated in pyrethroid resistant strains [41]. These suggest that the CYP9 family is induced by 
synthetic pyrethroids. In many cases, overexpression of P450s often leads to increased metabolism of 
pesticides, and increased enzyme activity often enhances detoxification [42]. In this experiment, 
both the molecular and biochemical data indicate that CYP9A61 is strongly induced by exposure to 
lambda-cyhalothrin, suggesting that CYP9A61 is a critical protein potentially associated with 
lambda-cyhalothrin metabolism. Further studies should be carried out to verify the function of 
CYP9A61 using RNAi to knock the transcript level down or metabolize lambda-cyhalothrin in vitro. 

In the present study, variations between CYP9A61 expression levels and PNOD activity post 
exposure were not precisely consistent with each other, especially for chlorpyrifos-ethyl. For 
chlorpyrifos-ethyl, the CYP9A61 transcript is 2.2-fold increased by up to 60 h exposure (Figure 4A), 
whereas PNOD activity was remarkably increased after 24 h of insecticide induction and weakly 
increased at 36 h, with subsequent activity repressed after 48 h treatment (Figure 6). Another peak 
appeared at 72 h. The observed differences between gene expression and enzymatic activity demonstrate 
that the inducibility of CYP9A61 transcripts is not precisely consistent with enzyme activity. A similar 
finding has been documented for CYP9Q1, CYP9Q2 and CYP9Q3 in A. mellifera [17]. We suspect that 
the reason for no direct correlation between PNOD enzymatic activity and single CYP gene expression 
could be that cytochrome P450-mediated (9-demethylation activity is under polygenic control. It's also 
worth pointing out that each species contains numerous P450s [43]. Although the number of 
P450s expressed in the codling moth is not certain, there are numerous P450s expressed in third-instar 
larvae. In the present work, only the expression level of CYP9A61 was investigated. In order to clarify 
the MFO activity variations correlated with CYP6A61 transcript variations, more P450 genes should be 
further investigated. This phenomenon could also be explained by Serikawa's viewpoint that there is not 
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always a correlation between transcript level changes and translation changes [44]. Gene expression is a 
complex process that is regulated at several levels and regulation can occur in the post-transcriptional 
step [45]. However, the increased CYP levels in mRNA transcripts and enzymatic activity may provide 
evidence in support of CYP9A61 potential involvement in the metabolism of chlorpyrifos-ethyl, and this 
may be a common detoxification mechanism for CYP9s induction in response to insecticides. The low 
increased ratios observed in both mRNA expression level and enzyme activity are due to the low 
induction dose selected. The differences of CYP9A61 response to these two insecticides may be 
explained by the differences in insecticides and by the treatment dosages used. 

4. Experimental Section 

4.1. Insects 

The apples with tunnels and frass produced in the skin after consumption by C. pomonella larvae 
were collected from abandoned apple orchards in Wuwei City (Gansu Province, China), a non-managed 
codling moth affected area. Low instar C. pomonella larvae were reared on the original collected whole 
apples. The higher instar larvae were removed from the apples and reared on flesh of the apple debris in 
the laboratory at 16:8 h light:dark photoperiod, 25 =t 1 °C and 60% ± 5% relative humidity [5] without 
exposure to any insecticides and chemicals. The larvae were reared until the adult stage for mating to 
obtain the next generation of neonates. This codling moth strain reared in the laboratory for 
6 generations was used in this study. 

4.2. Insecticides and C. pomonella Treatment 

Chlorpyrifos-ethyl and lambda-cyhalothrin with purities of >99% were obtained from Aladdin 
Reagent (Shanghai, China). These insecticides were dissolved in acetone. Rodriguez et al. demonstrated 
that the ZC50 of chlorpyrifos-ethyl and lambda-cyhalothrin on neonate codling moth were 157.65 mg/L 
and 0.35 mg/L [4], respectively, when the insecticides are applied on the surface of a 4 cm 
semi-artificial diet piece. In the present study, the topical application method [10] was used. 
The induction doses using topical application method for these two insecticides for codling moth have 
not been defined previously. A non-lethal dose, 12.5 mg/L for chlorpyrifos-ethyl and 0.19 mg/L for 
lambda-cyhalothrin was used as the inducing dose for the gene induced-expression and enzyme assay. 
Larvae were starved for 12 h, and a 1 |liL of chlorpyrifos-ethyl solution (12.5 ng) and lambda-cyhalothrin 
solution (0. 19 ng) was applied on the dorsum of each larva according to this topical application method [10] 
using an Eppendorf pipettor (Hamburg, Germany). Acetone was used for the control. After treatment, 
larvae were supplemented with the fiesh of apple pieces as food and kept under the conditions described 
above. Each experiment was performed in triplicate, and for each treatment in each time point, five 
third-instar larvae were used. Insecticide and acetone treated larvae were collected at 12, 24, 36, 48, 60 h 
for investigation of the mRNA expression level and enzymatic study. First to fifth instar larvae, male and 
female pupae and adults were sampled from untreated controls. Each assay was carried out on 3 groups of 
5 individuals. The larvae were not sexed due to this being unfeasible at that Ufe stage. Head, cuticle, silk gland, 
midgut, fat body, and Malpighian tubules were collected from 3 groups of 30 third instar larvae. Each 
group was flash fi*ozen in liquid nitrogen and kept at -80 °C for RNA extraction or enzyme preparation later. 
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4.3. RNA Extraction and cDNA Synthesis for Molecular Cloning 

The total RNA of C. pomonella was extracted from five third instar larvae using the RNAiso Plus Kit 
based on the manufacturer's instructions (Takara, Dalian, China). The quality and concentration of 
extracted RNA was examined by agarose gel electrophoresis and spectrophotometer analysis 
(Infinite M200 PRO, Tecan Group Ltd., Mannedorf, Switzerland). The RNA samples were then digested 
with DNase I (MB I) to remove the genomic DNA and for cDNA synthesis. For 3 '-rapid amplification of 
cDNA ends (3'RACE), first strand cDNA was synthesized from 1 |ig of total RNA using the SMART™ 
RACE cDNA Amplification Kit and MLV Reverse Transcriptase (Clontech, Dalian, China) as 
described by the manufacturer. For 5'-rapid amplification of cDNA ends (5'RACE), single-stranded 
cDNA was synthesized from 1 |Lig of total RNA by 5 '-Full RACE Kit (Takara, Dalian, China) according 
to the procedure recommended. The product was stored at -20 °C for future use. 

4.4. Molecular Cloning ofCYP9A61 cDNA by 3' and 5'RACE 

Degenerate primers used for 3' and 5 'RACE amplification were designed (Table 1) based on the 
conserved amino acid sequences of the insects' CYP9 family deposited in the NCBI database using the 
CODEHOP web tool [46] . 3'RACE PCR was conducted by nested PCR using P9-3F 1 and 1 Ox Universal 
Primer Mix (UPM, Clontech, Dalian, China) as the primer pair for the first round PCR reaction and 
P9-3F2 combined with lOx UPM for the second round amplification. For 5'RACE, PCR was carried out 
by nested PCR using P9-5R1 and Outer Primer as the primer pair for first round PCR amplification and 
followed by the second round amplification using P9-5R2 combined with Inner Primer (Table 1). 
The Outer and Inner Primers were supplied by 5'-Full RACE Kit. The PCR was conducted in a 
ClOOO Thermal Cycler (BioRad, Hercules, CA, USA). For both 3' and 5'RACE, the PCR conditions for 
the first round PCR reaction were: 94 °C for 3 min, followed by 35 cycles of 30 s at 94 °C, 30 s at 55 °C 
and 2 min at 72 °C, and then a final extension at 72 °C for 7 min. For the second round PCR, 
PCR reaction conditions were the same except the extension time was shortened to 1 min. The 
PCR products were separated on a 1% agarose gel electrophoresis and the bands with approximately 
1000 bp (for 3'RACE) and 700 bp (for 5'RACE) were gel purified using the Biospin Gel Extraction Kit 
(Bioer Technology Co., Ltd., Hangzhou, China) as described by the manufacturer. The purified 
fragments were ligated into the pMD-19 T vector (Takara, Dalian, China) and ten positive clones were 
sequenced (Shanghai Sunny Biotech Co., Ltd., Shanghai, China). Only one 3' and 5' cDNA end were 
obtained using those primers. The sequenced clones did not show any SNPs. Based on the obtained 
sequences of the 3' and 5'RACE, the putative full-length of CYP9A61 gene was amplified with specific 
primers P9F and P9R (Table 1). The PCR was conducted using high-fidelity Ex Tag polymerase 
(Takara, Dalian, China) to eliminate any potential error occurring by Tag DNA polymerase. 
The full-length cDNA sequence was subjected to the Open Reading Frame Finder (ORF Finder, 
http://www.ncbi.nlm.nih.gov/gorf/gorf.html) for open reading frame (ORF) prediction. The ORF of 
CYP9A61 cDNAs was amplified by using primer pairs P90RF-F and P90RF-R (Table 1). The PCR 
was performed at 94 ^C for 3 min, followed by 35 cycles of 94 ^C for 30 s, 55 ^C for 30 s and 72 ^C for 
2 min, with a final extension at 72 °C for 7 min. The PCR product was TA cloned and sequenced as 
described above. 
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4.5. Sequence Analysis 

The CYP9A61 amino acid sequence was deduced from ORF cDNA sequence using the ExPASy 
Proteomics web tool [47]. Using this web tool, the molecular weight (Mw) and theoretical isoelectric 
point ipl) was also predicted. The deduced amino acid sequence of CYP9A61 was aligned with 
Manduca sexta CYP9A5 (GeneBank number: AAD5 1038.1), H. virescens CYP9A1 (GeneBank 
number: AAC25787.1) and Bombyx mori (L.) CYP9A20 (GeneBank number: NP_00 1077079.1) 
which have been reported as possibly induced by insecticide exposure [21] using the 
ClustalW2 software [48]. The secondary structure prediction was performed with PSIPRED [49]. 
Motifs were predicted based on scanning the Prosite database [50] The presence of signal peptide was 
predicted by SignalP 3.0 [51] and transmembrane was determined using the Simple Modular 
Architecture Research Tool [52]. An un-rooted tree was constructed by MEGA 4 [53] under the 
Poisson-correction model using the neighbor-joining method with 1000 replications support. 

4.6. Real Time Quantitative PGR (qPCR) 

QPCR was carried out to study tissue distribution, stage specific and insecticide-induced expression 
of CYP9A61. The total RNA was extracted from collected tissues, the whole body of various 
developmental stages of C. pomonella, and insecticide and acetone (control) larvae. For developmental 
expression and insecticide induction expression assays, five individuals were used to isolate total RNA, 
using three independent extracts. For the tissue specific expression study, each tissue dissected from 
30 third instar larvae was regarded as an independent specimen to isolate total RNA, using three 
independent extracts. The first strand cDNA of each RNA sample for qPCR was synthesized using 
1 .5 |ig of total RNA using PrimeScripf^^ RT reagent Kit with gDNA Eraser (Takara, Dalian, China) in a 
20 |LiL reaction mixture as described by the manufacturer. The relative expression levels of CYP9A61 in 
different tissues, at each developmental stage, and under chlorpyrifos-ethyl and lambda-cyhalothrin 
exposure were assessed using the reference gene C. pomonella beta-actin (CpActin, GenBank number: 
KC832921) as a reference. In our previous study, CpActin was showed to be an excellent reference gene 
due to being expressed constantly in different tissues, developmental stages, and exposure to 
insecticides [54]. The primers used were QCYP9A61F and QCYP9A61R (Table 1) for target gene 
amplification, and QActinF and QActinR (Table 1) for endogenous control. The qPCR was performed 
on a BioRad iCycler iQ5 (BioRad, Hercules, CA, USA) using an UltraSYBR Mixture (CWBIO, Beijing, 
China) with 20 |liL reaction mixture consisting of 10 |liL UltraSYBR Mixture, 1.5 |iL cDNA, 0.1 |iL of 
each 10 |liM forward and reverse primers, and adding distilled water up to 20 |liL under the following 
cycling program: 1 cycle at 95 "^C for 10 min, 40 cycles at 95 °C for 15 s, 53 "^C for 1 min, and a melting 
curve program (55 to 95 °C in increments of 0.5 °C every 20 s) was conducted after PGR to confirm if 
the nonspecific product was produced. The standard curves of both CYP9A61 and CpActin were 
obtained using a linear gradient dilution (10-fold dilution) of cDNA as template. The amplification 
efficiencies, slopes, and correlation coefficients (R ) were 100.7%, -3.306 and 0.994 for CYP9A61 
standard curve, and 95%, -3.334 and 0.995 for CpActin, respectively. The relative CYP9A61 mRNA 
expression was calculated according to the 2~^^^^ method [55]. A negative control (NTC) reaction was 
carried out using RNase-free water instead of cDNA templates. A no transcription control (NRC) was 
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also conducted to eliminate DNA contamination. Each experiment contains three biological replicates, 
and each sample was technically repeated three times. 

4. 7. Enzyme Activity 

Whole bodies from insecticide-treated and acetone-treated third instar larvae were homogenized on 
ice for three replications of five individuals each in 2 mL ice-cold 0.1 M sodium phosphate buffer 
(pH 7.8) containing a final concentration of 1 mM phenylmethylsulfonylfluoride (PMSF), 1 mM ethylene 
diamine tetraacetic acid (EDTA) and ImM dithiothreitol (DTT). The homogenates were centrifuged at 
4 °C for 20 min at 12,000 g, and the supematants were used as an enzymatic source for the MFO assay [56]. 
The proteins were quantified by the Bradford method [57] using bovine serum albumin as the standard. 

The MFO activities were assessed by measuring PNOD activity on Nunc 96-well transparent 
microplates (Nunc, Roskilde, Denmark) as described by Yang et al, [56] with slight modification. 
The reaction mixture consisted of 180 |liL of 0.1 M sodium phosphate buffer (pH 7.8) containing 2 mM 
/?-nitroanisole, 10 |liL of enzyme. After 5 min of incubation at 30 °C, the reaction was initiated by adding 
10 |LiL of 9.6 mM reduced form of nicotinamide adenine dinucleotide phosphate (NADPH). A negative 
control reaction was conducted using heat inactivated enzyme (95 °C for 10 min) to eliminate any 
non-enzymatic activity. Moreover, a no-NADPH-added control was also performed. Pre-incubated 
sodium phosphate buffer was used for the blank. After 1 h of reaction at 30 °C, the absorbance was read 
in an Infinite M200 PRO Microplate Reader (Tecan, Mannedorf, Switzerland) at 405 nm. A standard 
curve was conducted with /?-nitrophenol, and the result was expressed as pmole /^-nitrophenol min"^ mg 
of protein~\ Each experiment contains three biological replicates and three replicates were conducted 
for each test and the average was used for statistical analysis. 

4.8. Statistical Analysis 

Statistical analysis was carried out by the SPSS statistics software (version 12, IBM, IBM Inc., 
Chicago, IL, USA). One-way analysis of variance (ANOVA, = 0.05) was used for statistical analysis 
to determine the significance of the difference in CYP9A61 mRNA expression for developmental stages 
and insecticide induction in the codling moth. The differences in CYP9A61 mRNA expression for tissue 
distribution was assessed by Student's Mest (*;? < 0.05; < 0.01; p < 0.001). Factorial ANOVA 
and Student's ^test were used to analyze the enzyme activity of C. pomonella treated with different 
insecticides and different durations of exposure to insecticide. The results were plotted using GraphPad 
Prism 5 (GraphPad Software, San Diego, CA, USA). 

5. Conclusions 

In summary, this is the first sequence of a P450 gene isolated from C. pomonella. We can conclude 
that CYP9A61 is a xenobiotic-metabolizing P450 gene. This study advances our understanding of 
CYP9A61 function, and regulation mechanism, and the CYP9A61 response to insecticide exposure. 
Further studies of CYP9A61 detoxification of xenobiotic compounds in vitro are required for accurate 
assessments of the metabolic functional activity toward these chemicals. 



Int. J. Mol Set 2013, 14 



24226 



Acknowledgments 

This research was supported by the Special Fund for Agro-scientific Research in the Public Interest of 
China (No. 200903042-03). We thank David Nelson for naming the P450 gene, and John Richard 
Schrock from Emporia State University (Emporia, KS, USA) for proofreading this manuscript before 
submission. We also thank Xi'en Chen, Ming Jiang and Long Ma for their help in the dissection of 
codling moth tissues. 

Conflicts of Interest 

The authors declare no conflict of interest. 
References 

1. Reuveny, H.; Cohen, E. Evaluation of mechanisms of azinphos-methyl resistance in the codling 
moth Cydia pomonella (L.). Arch. Insect Biochem. Physiol. 2004, 57, 92-100. 

2. Fuentes-Contreras, E.; Reyes, M.; Barros, W.; Sauphanor, B. Evaluation of azinphos-methyl 
resistance and activity of detoxifying enzymes in codling moth (Lepidoptera: Tortricidae) from 
central Chile. J. Econ. Entomol. 2007, 100, 551-556. 

3. Reyes, M.; CoUange, B.; Rault, M.; Casanelli, S.; Sauphanor, B. Combined detoxification 
mechanisms and target mutation fail to confer a high level of resistance to organophosphates in 
Cydia pomonella (L.) (Lepidoptera: Tortricidae). Pestic. Biochem. Phys. 2011, 99, 25-32. 

4. Rodriguez, M.A.; Marques, T.; Bosch, D.; Avilla, J. Assessment of insecticide resistance in eggs 
and neonate larvae of Cydia pomonella (Lepidoptera: Tortricidae). Pestic. Biochem. Phys. 2011, 
100, 151-159. 

5. Voudouris, C.C.; Sauphanor, B.; Frank, P.; Reyes, M.; Mamuris, Z.; Tsitsipis, J.A.; Vontas, J.; 
Margaritopoulos, J.T. Insecticide resistance status of the codling moth Cydia pomonella 
(Lepidoptera: Tortricidae) from Greece. Pestic. Biochem. Phys. 2011, 100, 229-238. 

6. Soleno, J.; Anguiano, O.L.; Cichon, L.B.; Garrido, S.A.; Montagna, CM. Geographic variability in 
response to azinphos-methyl in field-collected populations of Cydia pomonella (Lepidoptera: 
Tortricidae) from Argentina. Pest Manag. Sci. 2012, 68, 1451-1457. 

7. Sauphanor, B.; Bouvier, J.C.; Brosse, V. Spectrum of insecticide resistance in Cydia pomonella 
(Lepidoptera: Tortricidae) in Southeastern France. J. Econ. Entomol. 1998, 91, 1225-1231. 

8. Pasquier, D.; Charmillot, P.J. Effectiveness of twelve insecticides applied topically to diapausing 
larvae of the codling moth, Cydia pomonella L. Pest Manag. Sci. 2003, 60, 305-308. 

9. Reyes, M.; Frank, P.; Charmillot, P.J.; Loriatti, C; Olivares, J.; Pasqualini, E.; Sauphanor, B. 
Diversity of insecticide resistance mechanisms and spectrum in European populations of codling 
moth, Cydia pomonella. Pest Manag. Sci. 2007, 63, 890-902. 

10. Rodriguez, M.A.; Bosch, D.; Sauphanor, B.; Avilla, J. Susceptibility to organophosphate 
insecticides and activity of detoxifying enzymes in Spanish populations of Cydia pomonella 
(Lepidoptera: Tortricidae). J. Econ. Entomol. 2010, 103, 482^91. 



Int. J. Mol Set 2013, 14 



24227 



11. Sauphanor, B.; Cuany, A.; Bouvier, J.C.; Brosse, V.; Amichot, M.; Berge, J.B. Mechanism of 
resistance to deltamethrin in Cydia pomonella (L.) (Lepidoptera: Tortricidae). Pestic. Biochem. Phys. 
1997, 58, 109-117. 

12. Bouvier, J.C.; Bues, R.; Boivin, T.; Boudinhon, L.; Beslay, D.; Sauphanor, B. Deltamethrin 
resistance in the codling moth (Lepidoptera: Tortricidae): Inheritance and number of genes 
involved. Heredity 2001, 87, 456^62. 

13. Sanchez, D.M.; Wise, J.C.; Poppen, R.V.; Gut, L.J.; HoUingworth, R.M. Resistance of codling 
moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae), larvae in Michigan to insecticides with 
different modes of action and the impact on field residual activity. Pest Manag. Sci. 2008, 64, 
881-890. 

14. Scott, J.G. Cytochromes P450 and insecticide resistance. Insect Biochem. Mol. Biol. 1999, 29, 
757-777. 

15. Feyereisen, R. Insect Cytochrome P450. In Comprehensive Molecular Insect Science; Gilber, L.I., 
latrou, K., Gill, S., Eds.; Elsevier Science Publishers: London, UK, 2005; Volume 4, pp. 1-77. 

16. Li, X.C.; Schuler, M.A.; Berenbaum, M.R. Molecular mechanisms of metabolic resistance to 
synthetic and natural xenobiotics. Annu. Rev. Entomol. 2007, 52, 231-253. 

17. Mao, W.F.; Schuler, M.A.; Berenbaum, M.R. CYP9Q-mediated detoxification of acaricides in the 
honey bee {Apis mellifera). Proc. Natl. Acad. Sci. USA 2011, 108, 12657-12662. 

18. Schuler, M.A. P450s in plant-insect interactions. Biochim. Biophys. Acta 2011, 1814, 36^5. 

19. Zhou, J.L.; Zhang, G.R.; Zhou, Q. Molecular characterization of cytochrome P450 CYP6B47 
cDNAs and 5'-flanking sequence from Spodoptera litura (Lepidoptera: Noctuidae): Its response to 
lead stress. J. Insect Physiol. 2012, 58, 726-736. 

20. Nebert, D.W.; Gonzalez, F.J. P450 genes: Structure, evolution, and regulation. Annu. Rev. Biochem. 
1987, 56, 945-993. 

21. Rose, R.L.; Goh, D.; Thompson, D.M.; Verma, K.D.; Heckel, D.G.; Gahan, L.J.; Roe, R.M.; 
Hodgson, E. Cytochrome P450 (CYP)9A1 in Heliothis virescens: The first member of a new CYP 
family. Insect Biochem. Mol. Biol 1997, 27, 605-615. 

22. Poupardin, R.; Reynaud, S.; Strode, C; Ranson, H.; Vontas, J.; David, J.P. Cross-induction of 
detoxification genes by environmental xenobiotics and insecticides in the mosquito Aedes aegypti: 
Impact on larval tolerance to chemical insecticides. Insect Biochem. Mol. Biol. 2008, 38, 540-551. 

23. Bai, X.D.; Mamidala, P.; Rajarapu, S.P.; Jones, S.C.; Mittapalli, O. Transcriptomics of the Bed Bug 
(Cimex lectularius). PLoS One 2011, 6, el6336. 

24. Komagata, O.; Kasai, S.; Tomita, T. Overexpression of cytochrome P450 genes in 
pyrethroid-resistant Culex quinquefasciatus. Insect Biochem. Mol. Biol. 2010, 40, 146-152. 

25. Hasemann, C.A.; Kurumbail, R.G.; Boddupalli, S.S.; Peterson, J.A.; Deisenhofer, J. Structure and 
function of cytochromes P450: A comparative analysis of three crystal structures. Structure 1995, 
3,41-62. 

26. Zhou, X.J.; Ma, C.X.; Li, M.; Sheng, C.F.; Liu, H.X.; Qiu, X.H. CYP9A12 and CYP9A17 in the 
cotton boUworm, Helicoverpa armigera: Sequence similarity, expression profile and xenobiotic 
response. Pest Manag. Sci. 2010, 66, 65-73. 



Int. J. Mol Set 2013, 14 



24228 



27. Li, B.; Bisgaard, H.C.; Forbes, V.E. Identification and expression of two novel cytochrome P450 
genes, belonging to CYP4 and a new CYP331 family, in the polychaete Capitella capitata sp.I. 
Biochem. Biophys. Res. Commun. 2004, 525, 510-517. 

28. Stevens, J.L.; Snyder, M.J.; Koener, J.F.; Feyereisen, R. Inducible P450s of the CYP9 family from 
IbtvslI Manduca sexta midgut. Insect Biochem. Mol. Biol. 2000, 30, 559-568. 

29. Zhao, C.Q.; Tang, T.; Liu, J.Q.; Feng, X.Y.; Qiu, L.H. Identification and expression analysis of 
NADH-cytochrome bs reductase gene in the cotton boUworm, Helicoverpa armigera. Gene 2012, 
577, 96-102. 

30. Goff, G.L.; Hilliou, F.; Siegfried, B.D.; Boundy, S.; Wajnberg, E.; Sofer, L.; Audant, P.; 
French-Constant, R.H.; Feyereisen, R. Xenobiotic response in Drosophila melanogaster: 
Sex dependence of P450 and GST gene induction. Insect Biochem. Mol. Biol. 2006, 36, 674-682. 

31. Nikou, D.; Ranson, H.; Hemingway, J. An adult-specific CYP6 P450 gene is overexpressed in a 
pyrethroid-resistant strain of the malaria vector. Anopheles gambiae. Gene 2003, 575, 91-102. 

32. Snyder, M.J.; Stevens, J.L.; Andersen, J.F.; Feyereisen, R. Expression of cytochrome P450 genes 
of the CYP4 family in midgut and fat body of the tobacco hornworm, Manduca sexta. 
Arch. Biochem. Biophys. 1995, 321, 13-20. 

33. Arrese, E.L.; Soulages, J.L. Insect fat body: Energy, metabolism and regulation. Annu. Rev. Entomol. 
2010, 55, 207-225. 

34. Chung, H.; Sztal, T.; Pasricha, S.; Sridhar, M.; Batterham, P.; Dabom, P.J. Characterization of 
Drosophila melanogaster cytochrome P450 genes. Proc. Natl. Acad. Sci. USA 2009, 106, 573 1-5736. 

35. Li, X.C.; Schuler, M.A.; Berenbaum, M.R. Jasmonate and salicylate induce expression of herbivore 
cytochrome P450 genes. Nature 2002, 419, 712-715. 

36. Liu, N.N.; Li, T.; Reid, W.R.; Yang, T.; Zhang, L. Multiple cytochrome P450 genes: 
Their constitutive overexpression and permethrin induction in insecticide resistant mosquitoes, 
Culex quinquefasciatus. PLoS One 2011, 6, e23403. 

37. Bouvier, J.C.; Boivin, T.; Beslay, D.; Sauphanor, B. Age- dependent response to insecticides and 
enzymatic variation in susceptible and resistant codling moth larvae. Arch. Insect Biochem. 2002, 
57, 55-66. 

38. Harrison, T.L.; Zangerl, A.R.; Schuler, M.A.; Berenbaum, M.R. Developmental variation in 
cytochrome P450 expression in Papilio polyxenes in response to xanthotoxin, a hostplant 
allelochemical. v4rc/z. Insect Biochem. 2011, 48, 179-189. 

39. Strode, C; Wondji, C.S.; David, J.P.; Hawkes, N.J.; Lumjuan, N.; Nelson, D.R.; Drane, D.R.; 
Karunaratne, S.H.P.P.; Hemingway, J.; Black, W.C., I.V.; et al. Genomic analysis of detoxification 
genes in the mosquito Aedes aegypti. Insect Biochem. Mol. Biol. 2008, 38, 1 13-123. 

40. Gonzalez, I.F.; Quinones, M.L.; Lenhart, A.; Brogdon, W.G. Insecticide resistance status of 
Aedes aegypti (L.) from Colombia. Pest Manag. Sci. 2011, 67, 430-437. 

41 . Yang, Y.; Yue, L.; Chen, S.; Wu, Y. Functional expression of Helicoverpa armigera CYP9A12 and 
CYP9A14 in Saccharomyces cerevisiae. Pestic. Biochem. Phys. 2008, 92, 101-105. 

42. Denison, M.S.; Whitlock, J.P. Xenobiotic-inducible transcription of cytochrome P450 genes. 
J. Biol. Chem. 1995, 270, 18175-18178. 

43. Scott, J.G. Insect Cytochrome P450s: Thinking beyond Detoxification. In Toxicology and 
Molecular Biology', Liu, N.N., Ed.; Research Signpost: Kerala, India; 2008; pp. 1 17-124. 



Int. J. Mol Set 2013, 14 



24229 



44. Serikawa, K.A.; Xu, X.L.; Mackay, V.L.; Law, G.L.; Zong, Q.; Zhao, L.P.; Bumgamer, R.; 
Morris, D.R. The transcriptome and its translation during recovery from cell cycle arrest in 
Saccharomyces cerevisiae. Mol Cell. Proteomics 2003, 2, 191-204. 

45. Mata, J.; Marguerat, S.; Bahler, J. Post-transcriptional control of gene expression: A genome-wide 
perspective. Trends Biochem. Sci. 2005, 30, 506-514. 

46. Rose, T.M.; Henikoff, J.G.; Henikoff, S. CODEHOP (Consensus-Degenerate Hybrid 
Oligonucleotide Primer) PCR primer design. Nucleic Acids Res. 2003, 31, 3763-3766. 

47. Gasteiger, E.; Gattiker, A.; Hoogland, C; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: 
The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 
3784-3788. 

48. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; 
Valentin, F.; Wallance, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. 
Bioinformatics 2007, 25, 2947-2948. 

49. McGuffin, L.J.; Bryson, K.; Jones, D.T. The PSIPRED protein structure prediction server. 
Bioinformatics 2000, 16, 404-405. 

50. Hulo, N.; Bairoch, A.; Bulliard, V.; Cerutti, L.; Castro, E.D.; Langendijk-Genevaux, P.S.; Pagni, M.; 
Sigrist, C.J.A. The PROSITE database. Nucleic Acids Res. 2006, 34, D227-D230. 

51. Bendtsen, J.D.; Nielsen, H.; von Heijne, G.; Brunak, S. Improved prediction of signal peptides: 
SignalP 3.0. J. Mol. Biol. 2004, 340, 783-795. 

52. Schultz, J.; Milpetz, F.; Bork, P.; Pouting, CP. SMART, a simple modular architecture research 
tool: Identification of signaling domains. Proc. Natl. Acad. Sci. USA 1997, 95, 5857-5864. 

53. Tamura, K.; Dudley, J.; Nei, M.; Kumar, S. MEGA4: Molecular evolutionary genetics analysis 
(MEGA) software version 4.0. Mol. Biol. Evol. 2007, 24, 1596-1599. 

54. Yang, X.Q.; Zhang, Y.L. Molecular cloning, sequence analysis and mRNA expression stability 
detection of p-actin gene in codling moth Cydia pomonella (L.). J. Northwest A&F Univ.: 
Nat. Sci. Ed. 2013, in press. 

55. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative 
PCR and the 2""^^^^ method. Methods 2001, 25, 402^08. 

56. Yang, Y.H.; Wu, Y.D.; Chen, S.; Devine, G.J.; Denholmb, L; Jewess, P.; Moores, G.D. 
The involvement of microsomal oxidases in pyrethroid resistance in Helicoverpa armigera from 
Asia. Insect Biochem. Mol. Biol. 2004, 34, 763-773. 

57. Bradford, M.M. A rapid and sensitive method for the quantification of microgram quantities of 
protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248-254. 



© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



